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Abstract 

The study of exclusive B decays in perturbative QCD are complicated by the end- 
point problem. In order to perform the perturbative calculation, the Sudakov effects 
are introduced to regulate the endpoint singularity. We provide a systematic anal- 
ysis with leading and next-to-leading twist corrections for B — > n form factors in 
pQCD approach. The intrinsic transverse momentum dependence of hadronic wave 
function and threshold resummation effects are included in pQCD approach. There 
are two leading twist B meson distribution amplitudes (or generally wave functions) 
in general. The QCD equations of motion provide important constraints on B me- 
son wave functions. The reliability of pQCD approach in B — * n form factors is 
discussed. 70% of the result comes from the region a s (t)/ir < 0.2 and 38% comes 
from the region where the momentum transfer t > lGeV. The conceptual problems 
of pQCD approach are discussed in brief. Our conclusion is that pQCD approach 
in the present form cannot provide a precise prediction for B — > it transition form 
factors. 
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1 Introduction 



The exclusive B decays provide an important test of the standard model of particle physics. 
With the running of BaBar and Belle B-factories and the proceeding of future B physics 
projects ( BTeV and LHCb etc.), a large amount of B mesons will be accumulated to 
explore the origin of CP violation and determine the CKM parameters, such as the angles 
a, j3 and 7 in unitary triangle. In most cases, the limitation of our theoretical ability 
prevents the precise prediction for exclusive B decays so we have to refer to some phe- 
nomenological approaches. 

The large B meson mass m B establishes a large scale so that perturbative QCD 
(pQCD) may be applicable in exclusive B decays. Recently, two different approaches 
based on perturbative QCD were proposed to calculate the exclusive B decays. One 
approach, usually say QCD factorization approach, states that two body hadronic B 
decays can factorize and the relevant amplitude can be written as the convolutions of 
non-perturbative quantities (the light-cone distribution amplitudes of mesons and semi- 
leptonic form factors) and perturbatively calculable hard scattering kernels [1]. The other 
approach is the modified pQCD, or say pQCD approach for simplicity. In pQCD ap- 
proach, the semi-leptonic form factors for B — > P transition (where P represents light 
meson) are claimed to be perturbatively calculable [2]. 

The central difference between these two approaches is whether B — > P transition 
form factors are perturbatively calculable, or specifically, whether Sudakov effects can 
cure endpoint singularity. In B — > P transition form factors, the endpoint singularity 
is generated even at leading twist. In pQCD approach, the transverse momentum are 
retained to regulate the endpoint singularity and Sudakov double logarithm corrections are 
included to suppress the long distance contributions from configuration of large transverse 
separation. In [3] , the authors investigated the reliability of Sudakov effects in B — > ir form 
factors. Their conclusion is that Sudakov suppression are so weak that they can not be 
applied in B decays. If their criticism is right, it would bring the entire thought about 
Sudakov effects into crisis. 

In our previous study [4], we investigated the Sudakov effects in QCD factorization 
approach. Our conclusion is contrary to the above criticism. Sudakov effects play an 
important role in exclusive B decays. It is well known that the standard approach or 
say BLER approach [5] is questionable at experimentally accessible energy scales, typi- 
cally a few GeV region [6]. This question also occurs in the hard spectator scattering in 
B — > P\P 2 decays. The QCD factorization approach gives a finite result at leading twist. 
Carefully study shows that there are large contributions coming from the region where 
the momentum transfers are small. The introduction of Sudakov suppression enlarges 
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the application range of pQCD and makes it self-consistent [7]. Although the relevant 
numerical result is power suppressed, the Sudakov effects is potentially important. We 
disagree with the comment that the main purpose of introduction of Sudakov suppression 
is to improve the computation accuracy [3]. At twist-3 level, the chirally enhanced power 
correction is logarithmically divergent, so the QCD factorization approach can not apply. 
A phenomenological parametrization method is designed to regulate the endpoint singu- 
larity [8]. This method is not self-consistent and it introduces an arbitrary infrared cutoff. 
Sudakov suppression establish a natural infrared cutoff. Obviously, pQCD approach has 
the advantage of the parametrization method. 

The disagreement of conclusions in [2, 4] and [3] motivates us to reconsider the relia- 
bility of pQCD approach in B — > n form factors. In this paper, we present a systematic 
study of B — > 7r form factors in pQCD approach and examine the reliability of pQCD 
calculation. Compare with the recent pQCD analysis of B — > n form factors [9], the 
present analysis contains three new theoretical ingredients: 

• The intrinsic transverse momentum dependence of wave functions for B and pion 
mesons are included. The importance of intrinsic transverse momentum dependence 
is first pointed out in pion form factor in [10]. Because Sudakov suppression is not 
severely strong for real mg, the effect of intrinsic transverse momentum dependence 
of wave function can not be neglected. 

• B meson contains two leading wave functions. The assumption of single B meson 
wave function in the previous pQCD analysis is not valid. Equations of motion in 
HQET provide important constraint on the choice of B meson wave functions. 

• The threshold resummation effect for B meson is taken into account. The perturba- 
tive analysis depends on the endpoint behavior of B meson distribution amplitudes. 
The jet function obtained from threshold resummation suppress the endpoint con- 
tribution, so it modifies the power behavior of endpoint contribution. We expect 
that it can improve the perturbative result since one of the two B meson distribution 
amplitudes does not vanish at endpoint. 

In addition, the chirally-enhanced power corrections (twist-3 contribution) are included 
in our analysis. It is noted that the contribution of chirally-enhanced power corrections 
is comparable with or even larger than leading twist contribution [9]. Our study confirm 
it. The numerical results of B — > 7r form factors at large recoil region in pQCD approach 
are consistent with derivations of QCD sum rules [11]. 

The remainder of the paper is organized as follows. Sec. 2 introduces the theoretical 
ingredients of pQCD approach. In sec. 3, we present the formulas for B — > n form factors 
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in pQCD approach, perform the numerical analysis, and examine the reliability of pQCD 
approach. Finally, in sec. 4, the conclusions and discussions are presented. 

2 The B — > it form factors in pQCD approach 

B — > 7i form factors are the basic non-perturbative parameters in semi-leptonic B — > ixlv 
decays and exclusive, nonleptonic decays such as B — > 7T7t, 7T.K" decays in QCD factoriza- 
tion approach. While in pQCD approach, it can be perturbatively calculated from the 
universal hadronic distribution amplitudes or wave functions. In this section, we will give 
a general discussion of B — > ix form factors and introduce the main ingredients of pQCD 
approach. 

The form factors of B — > n are defined by the following Lorentz decompositions of 
biquark current matrix elements: 

( 7 r(P 7r )|fZ 7M 6|S(P B )) = (P B + P n - ^^q),F^(q 2 ) + ^^^(g 2 ), (1) 

where q = Pb — P-k is the momentum carried by lepton pair in semi-leptonic B — > 7r/z/ 
decays. At large recoil limit, q 2 = 0, the two form factors F+* , F Bn reduce to one 
parameter, F Bn {0) = F Bn {0). 

First, we give our conventions on kinematics. We work in the rest frame of B meson. 
The mass difference of b quark and B meson is negligible in the heavy quark limit and 
we take approximation m h = m B in our calculation. The masses of light quarks u, d 
and 7i meson are also neglected. For discussion, the momentum is described in terms of 
light-cone variables. We take k = (^, k±) with /c ± = k° ± k 3 and k± = (k 1 , k 2 ). The 
scalar product of two arbitrary vectors A and B is A • B — A^B^ = A+B + A B+ —A±-B±. 
The momentum of pion is chosen to be in the " — " direction. Under these conventions, 

p B = ( 7f ' P * = (°' p " = (^Vt' °' ^ with ^ = 1 - * = & B We define 

two light-like vectors n + = (\/2, 0, 0_l) and n + = (0, v^2, 0j_). 

We denote £ as the momentum fraction of spectator anti-quark in B meson, and x as 

the momentum fraction of the anti-quark in pion. As plotted in Fig .1. 

2.1 The physical picture of factorization in B — > n form factors 

The most important ingredient of pQCD is factorization, i.e., the separation of the long- 
distance dynamics from the short distance dynamics which is pertubatively calculable. 
Although rigorous proof of the factorization theorem is technically intricate, the physical 
picture is simple and intuitive. 



4 



First, we illustrate the factorization of pion electromagnetic form factor in order to 
introduce the basic idea of standard approach and the modified approach. A highly 
virtual photon collides on a quark in initial pion and makes it change its direction. In 
exclusive process 7*7r — > ir, the simple power counting shows that at large momentum 
transfers, the valence quarks dominate the process and the higher Fock states and the 
intrinsic transverse momentum is power suppressed. The large momentum Q means the 
high resolution at the distance of 1/Q. In a small distance region, a parton only sees the 
other parton relative to it and the initial and final pion can be considered as a small- 
size color-singlet object during the hard interaction. What is important for small-size 
color-singlet pion is that soft gluon corrections cancel at the leading twist. This is the 
well-known phenomena, "color transparency" . The long-distance interactions only occur 
before and after the hard interaction, so their effects can be factorized into the initial and 
final pions respectively. Since the hard interaction is restricted in the short distance, we 
only need the probability for the qq pair in pion to be within the transverse distance of 
l/Q, in other words, the distribution amplitude of pion (f)(x,Q 2 ). The scale Q acts as 
factorization scale as well as renormalization scale. The above discussions can be grouped 
into the standard factorization formula for electromagnetic form factor of pion 



The standard factorization is proved successfully in the asymptotic limit. The appli- 
cation of this standard approach at the experimentally accessible energy scales, typically 
a few GeV region, is criticized in [6]. The authors pointed out that the perturbative 
calculations for electromagnetic form factor have a large amount of contributions coming 
from soft endpoint region (x, y — > 0) where the perturbative analysis is invalid. This is 
the well-known "endpoint" problem even there is no endpoint singularity in the convo- 
lution of total amplitude. The Sudakov effects are introduced to modify the endpoint 
behavior and make pQCD applicable in few GeV region [7]. The basic idea of this modi- 
fied approach, or say pQCD approach is using the mechanism of Sudakov suppression to 
suppress long-distance contributions of large transverse separations. Sudakov suppression 
establishes a factorization scale 1/6 in addition to the momentum scale Q. The pion with 
small transverse separation has small color dipole moment so that QCD factorization is 
revised. In pQCD approach, Sudakov suppression plays a crucial role. 

Second, We discuss the factorization in B — > it form factors in the large uib limit. 
B meson is a heavy-light system. It has large size in longitudinal as well as transverse 
directions. If the soft B meson wave functions overlap with the wave functions of final 
meson, the separation of long-distance and short distance dynamics is impossible. From 
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this point, we conclude that B — > Z)W from factors in heavy quark limit are soft dominant 
and cannot be calculated by pQCD. For B — > 7r form factors, the situation is different. 
The pion carries the energy of ^ when q 2 = 0. When the fractional momentum of the 
antiquark in pion is far away from the endpoint region, i.e. x ^> 0, the pion is highly 
restricted in longitudinal direction because of Lorentz contraction. The soft gluons that 
attach to small size color-singlet pion decouple in the large m# limit. So factorization is 
applicable. Simply to say, the soft spectator anti-quark in B meson must undergo hard 
strong interaction to change into the fast moving parton in pion. 

The above perturbative picture is destroyed by the endpoint contribution where anti- 
quark in pion carries the momentum of order Aqcd- In the standard approach, B — > 
7r form factor is logarithmically divergent at leading twist (twist-2 for our case). The 
divergence becomes more serious for next-to-leading twist (twist-3) contribution and the 
linear divergence occurs. The origin of this divergence is that the distribution amplitude 
does not provide enough suppression at endpoint. In order to apply the perturbative 
picture, the soft endpoint contribution must be suppressed. As we have discussed for 
pion electromagnetic form factor, Sudakov effects can provide such suppression. With the 
help of Sudakov suppression, pQCD factorization can be applicable. In pQCD approach, 
the B — > 7r form factors can be generally written as convolutions of wave functions and 
hard scattering kernels 

F+,o = Jd^dx d 2 b B d% B (f , 6 B , b n , n) T + ,o(£, x, b B , b v , /i). (3) 

2.2 Pion wave functions 

In the standard approach, the transverse momentum is assumed power suppressed and 
neglected. When the two valence quarks in pion carry large longitudinal momentum, 
i.e., the momentum fractions x, x ^> 0, the transverse momentum is power suppressed 
compared to longitudinal momentum. In the endpoint region, the transverse momentum 
is important and cannot be neglected. In pQCD approach, transverse momentum is 
retained in hard scattering kernels and the non-perturbative parameters are hadronic 
wave functions in general. 

The two valence quarks in pion have transverse momentum as well as the longitudinal 
momentum. Compared to the collinear limit, the momentum of the quarks in pion ( with 
momentum P n ) changes to 

k\ = xpn + k±, k 2 = xp n - k±, (4) 

where x and x denote the longitudinal momentum fractions of quark and anti-quark, 
respectively. For our case, the meson is on-shell and the partons are slightly off-shell. 
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The off-shellness of the parton is proportional to k\ which is power suppressed compared 
to m 2 B . 

As we have discussed above, the intrinsic transverse momentum dependence of wave 
functions of pion and B meson can not be neglected because Sudakov suppression is not 
strong in the real tub energy. The distribution amplitude can be obtained from the 
integration of wave function over the transverse momentum 

<f>(x)= I d 2 k ± ^(x,k ± ), (5) 

J\k ± \<ix 

where /i is the ultraviolet cutoff. 

Similar to the definition of distribution amplitudes with leading and next-to-leading 
twist [12], the light-cone wave functions of pion are defined in terms of bilocal operator 
matrix element 

(n(p)\qp(z)q a \0) = ^ £ dx J d 2 k^ z ~ % ^ (6) 



X { $l5^n(x, k ± ) ~ V p (x, k±) - a^z" 
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/ ) a/3 

where /„- is the decay constant of pion. The parameter /j, n = m 2 /(m u + m^) for charged 
pion. For neutral pion, we use the same parameter as the charged pion. ty n , ^ p and 
ty a are the twist-2 and twist-3 wave functions, respectively. The twist-3 wave functions 
contributes power corrections. But at energy scale, the chirally enhanced parameter 
r\ = ^ ~ is n °t small. So the twist-3 contribution should be considered in B 

decays. 

In order to perform the analysis in pQCD approach, we need to transform the param- 
eters in terms of coordinate variable in Eq.(6) into the momentum space configuration. 
We use the momentum projection given in [13]: 

where % = . The wave functions fy^, ^p, may have different transverse 

momentum dependence, this will make the calculation difficult. In order to simplify dis- 
cussions, we assume the same transverse momentum dependence for these wave functions. 

In pQCD approach, the convolutions of wave functions and hard scattering kernel are 
presented in transverse configuration 6-space. We need to define wave function in 6-space 
through Fourier transformation 



(x,b) = J d 2 k ± e~ ii;± -h(x,k ± ). (8) 
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From Eq.(5) and (8), we can obtain 

<f>(x)=V(x,b = 0). (9) 

The impact parameter b represents the transverse separation of quark and anti-quark in 
pion. Eq.(9) shows that the distribution amplitude is equal to the wave function at zero 
transverse separation. 

The intrinsic transverse momentum dependence of wave function of pion is unknown 
from the first principle in QCD. We take a simple model in which the dependence of the 
wave function on the longitudinal and transverse momentum can be separated into two 
parts: 

V(x,k ± ) = (f>(x) x £(Jfej_), (10) 

where <f>(x) is the pion distribution amplitude. <f>(x) and S(fcj_) satisfy the normalization 
conditions 

j\(x) = l, Jd 2 k ± Z(k ± ) = 1. (11) 

The k± dependence of the wave function is contained in H(k±). In [10], T,(k±) is assumed 
to be a Gaussian distribution, 
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E(fc±) = ^eq>(-/Ffcl), (12) 

71 

Transforming it into the transverse configuration b-space, we obatin 

4(3 

The oscillating parameter (3 is fixed by requiring the root mean square transverse momen 
turn (r.m.s.), (k^) 1 / 2 is the order of Aq C d- Their relation can be obtained from 



E(6) = / d 2 k ± e -^- 5 E(fc ± ) = exp(-— ), (13) 



, tidxfd 2 k ± kl\^(x,k A ^ 2 



X / 1 da;/tf 2 jfej.|*(a;,A;j.)| 2 2/3 2 ' 

Thus, (k 2 } 1 / 2 = -j^. If the root mean square transverse momentum {k 2 ^ 1 / 2 = 0.35GeV, 
then (3 2 = 4GeV~ 2 . The detailed discussions about choice of the oscillation parameter (3 
can be found in [10]. 

2.3 B meson wave functions 

The intrinsic dynamic of B meson is different from the case of light pion meson. The 
momentum components of the spectator quark I are of order Aqcd- The most convenient 
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tool to describe B meson is heavy quark effective theory (HQET). Since we have chosen 
P n in the " — " direction, the hard scattering amplitude does not dependent on /_, the 
/_ dependence of the wave functions can be integrated out. In HQET, the B meson 
wave functions are defined by the general Lorentz decomposition of the light-cone matrix 
element [14, 3] 



(0\qp(z)b a (0)\B(p B )) 



ifB \i> B J rm B 



2^+(z 2 ,t) + BK - ' By - ' f 



t 



75 



,(15) 



a/3 



with v = ^ and t = v ■ z. A path-ordered exponential is implicitly present in the 
gauge-independent matrix element. 

The wave functions in momentum space can be obtained through the Fourier trans- 
formation: 



(16) 



where / = (l + /y/2, 0, l±) and z = (0, z_/-\/2, z±), so that z 2 = —z]_ and t = z_/2. 

A momentum projection for the matrix element of B meson is needed to simplify the 
calculation. In [13], the authors obtain a projection which is valid for distribution ampli- 
tudes. We extend their analysis to include the transverse momentum dependence of wave 
functions. We observe that their formulas need no modification except the replacement of 
distribution amplitudes by wave functions. The proof of this point is given in Appendix. 
So, the momentum-space projection operator for B meson is: 



M B _ ifB \ j>B + m B 
Map ~ -T \ ~ 



^ + + *B*- ~ A(J+, W-w 



d_ 



75 



(17) 



a/3 



where A(Z + , l ± ) = jf dl(^ B (l, l ± ) - ¥+(J, l ± )). 

The projection operator can be represented in the form which is helpful to compare 
with the results in the previous pQCD analysis [9] 



M a/3 



ift 



(&B + m B) 



75 



(18) 



a(3 



where ^b, ^ b are defined by 



2 1 2 
In [1], the authors consider a different prjection operator for B meson 



{{i> B + m B ) lb [^Bi + i^B2]} , 



(19) 



(20) 
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where ^bi,^B2 are defined by 

*B1 = *£, ^£2 = ^~ V (21) 

The above three projection operators are equivalent up to leading power in 1/m^. The 
projection given in Eq.(20) neglects the A term which is proportional to l + /m b and thus 
power suppressed. 

The B meson distribution amplitudes are obtained from transverse momentum integral 
of the relevant wave functions or from wave functions at zero transverse separation. 0^ 
and 4> B ,4> B are distribution amplitudes relative to wave functions of and ^I B ,^/ B 
respectively and they are related by 

<PB = ~ , <?B = (22) 

The equations of motion impose constraint on the wave functions. Using the equation 
of motion for the light quark of B meson and neglecting the effects of three-parton and 
higher Fock states, we can obtain [13, 3] 



dt t 



dz 2 4 dt 2 



= 0, (23) 

z 2 =0 

= . (24) 

z 2 =0 



Similar to the discussion for pion wave function, we consider a model in which the 
dependence on the longitudinal and transverse momenta of B meson is factorized: 



*b(&J±)=^(0xEb(J±), (25) 
where £ = ^ is the longitudinal momentum fraction. The 0^ are the two distribution 



amplitudes of B mesons. The normalization conditions are 

J d£0|(O = 1, J d 2 t ± J: B (l ± ) = 1. (26) 

Under the above assumptions, the two constraints become 



<PU0 = -^(0, (27) 
£ 2 ™!0b(O = "1^(0, (28) 

where u% = J <Pl±Pjl%(l±). 



10 



The distribution amplitudes of <p B can be solved analytically and the corresponding 
solution is: 




'2fm| , em] 

4 

The parameter uj b is fixed by the value of root mean square transverse momentum 
of B meson. For the transverse momentum dependent function, we choose the Gaussian 
distribution: 

The parameters ojb and uj' b are not independent, they are related by ujb = V^co'b- The 
root mean square transverse momentum (l^) 1 ^ 2 = ^| = uj' b . In transverse configuration 
b-space, the transverse momentum dependent function is 

£+(M = exp(-^M), (31) 

where bs is the conjugated variable of l±. For £^(&s) function, we assume it is the same 
as T,%(b B ). 

A more rigorous study of B meson wave functions in HQET is given recently in [15]. 
The authors find an analytic solution in the heavy quark limit using the QCD equations of 
motion. The transverse momentum distribution can also be derived without the assump- 
tion of Gauss behavior. The obtained B meson wave functions is model-independent. 

2.4 Sudakov form factor 

There are two types of resummation: Sudakov resummation (or say 6-space resumma- 
tion) and threshold resummation. These two resummation effects lead to suppression in 
different space: the region with large transverse separations b for Sudakov resummation 
and the small longitudinal fractional momentum x region for threshold resummation. 

First, we discuss the Sudakov resummation in brief. At a s order, the overlap of soft 
and collinear divergences produce double logarithms — c\v? Qb. The transverse impact 
parameter b is used to regulate the infrared divergence. In transverse configuration b space, 
the Sudakov double logarithms are resumed up to next-to-leading-log approximation. A 
exponential factor e~ s ^ x,b ^ will be obtained from 6-space resummation. we present the 
explicit expression of the exponent s(x, b, Q) appearing in Sudakov form factor. 

Define the variables, 

q = In (xQ/(V2A QCD )) , b = \n(l/bA QCD ). (32) 
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According to [16], the exponent s(x,b,Q) is presented up to next-to-leading-log approxi- 
mation 



s(x, b, Q) = 
(q 

q In - 

2/3i U 



A(1) (~ t\ A{2) {q A 



H — 



4/3? 
AW fa 



2(3 1 V / 4/3? 
ln(2g) + 1 ln(26) + 1 



4/3? 4A 



In 



In 



+ 



In 



3 2 7B -1' 



8/3? ""'V 2 
AP)/3 2 g - 6 



ln(2g) 



8/3? 
ln(26) + 1 



ln 2 (2g) - ln 2 (26) 



16/3? 62 



2 ln(26) + 1 



q b 
A^ (31 q - b 



AW(3 2 
16/3? 



21n(2g) + 3 21n(26)+3 



+ 



432/3? 6 3 



91n 2 (26) + 61n(26) + 2 



AWft 
1728/3? 



18 ln 2 (2g) + 30 ln(2g) + 19 18 ln 2 (26) + 30 ln(26) + 19 



b 2 



(33) 



where the coefficients /3j and A^ are 

33 - 2n f 



Pi 



12 



02 



3 



153 - 19n f 
24 
- 2 10 



IE 



A( 2 ) = ^----n / + -/3 1 ln( — 
9 3 27 ; 3 V 2 



(34) 



with the Euler constant. 

The running coupling constant a s up to next-to-leading-log is written as 

1 /3 2 mln0it7A 2 ) 



7T ~ /3i ln(/x 2 /A 2 ) /3? InV/A 2 ) 



(35) 



The exponent s(x,b,Q) is obtained under the condition that xQ/y/2 > 1/6, i.e. the 
longitudinal momentum should be larger than the transverse degree. So s(x, b, Q) is 
defined for q > b, and set to zero for q < b. The previous formulas [17] about the 
exponent s(x, b, Q) picks up the most important first six terms in the first and second 
lines of the expression of s. Note that the sign of the the fifth and sixth terms are different 
from those in [17]. 

The Sudakov form factor factor e~ s ( x ' b '^ falls off quickly in large b region and vanishes 
as b > 1/Aqcd- Therefore it suppresses the long-distance contribution, which is called 
Sudakov suppression. In axial-gauge, the Sudakov form factor is included in each hadronic 
wave function. So we can define the pion and B hadronic wave functions with Sudakov 
corrections as 



K- 1) = exp(-S , ff )*°(x, b v , t), tf B (f , b B ; t) = exp(-S B )V° B (£, b B , t), (36) 
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where t is the factorization scale in the hard scattering kernel. The ty®(x, b n , t) and 
^/° B (^,bB,t) are wave functions without Sudakov corrections. 

Combining with the evolution of wave functions and hard scattering kernel, a complete 
factor e~ S7T < B for pion and B mesons can be given as 



About the Sudakov form factor e s ( x > b < m B) f or b and pion meson, some comments are 
in order: 

• The Sudakov factor for B meson is only associated with the light quark since there 
is no collinear divergence associated with the heavy b quark. Due to the fact that 
the momentum of light spectator quark is concentrated at the order of Aq CD , we 
may expect the Sudakov effects is small because of the suppression of B meson wave 
function. Our numerical analysis shows that its effect is at one percent level. 

• Notice that the effect of the evolution term in the last of eq.(37) gives slight en- 
hancement. If Sudakov suppression takes place, the enhancement will be masked ( 
for example, the case of pion ). On the other hand, if Sudakov suppression is not 
effective ( for example, the case of B meson), this enhancement will emerge. 

• In the endpoint region, the exchanged gluon carries small longitudinal momentum. 
The transverse degree cannot be neglected. The mechanism of Sudakov suppression 
begin to take into effect, and the long-distance dynamics from large transverse 
separations are suppressed by Sudakov effects. If the momentum fraction of one 
quark x is small, the Sudakov form factor associated with it is 1, no suppression 
takes place. However the Sudakov form factor associated with the other quark will 
provide strong suppression because the momentum fraction of this quark x is large 
now. 

2.5 Threshold resummation 

Now, we discuss the threshold resummation effects. Double logarithms a s In 2 x can be 
produced by the higher order loop corrections. It diverges at the endpoint. If the end- 
point contribution is important, the double logarithms a s \v? x needs to be resumed to 
all orders. We review the basic idea of threshold resummation. Our discussion underlies 
on the analysis given in [9, 18]. The vertex corrections at a s order produce the double 




1 ln(t/A QCD ) 
Pi ln(l/(^A QCD ))' 



(37) 
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logarithms — ^-Cp In 2 x where Cp = 4/3 is the color factor. This collinear divergences can 
be factorized into a quark jet function S t (x). In order to resume the double logarithms 
to all orders, it is necessary to introduce the moment (N) space. In N space, the Sudakov 
factor has the exponential form up to the accuracy of leading-log approximation (LL), 

^ (LL) (iV)=exp[-i 7 ^ ) ln 2 iv], (38) 

where the anomalous dimension = ol s Cf/t^- 

The jet funation S t (x) can be obtained from S t (N) through the transfromation 

ra+ioo A AT , . 

S t (x) = f-(l - x)- N S t (N)S?\N) , (39) 

J a— ioc Z7TI 

where a is an arbitrary real constant larger than all the real parts of poles involved in 
the integrand. The sI°\n) comes from Mellin transformation of the initial condition 

Si°\x) = l, 

S?\N) = jT 1 dx(l - xf-'S^ix) = 1 . (40) 

The upper index "0" means that there is no QCD corrections. 
The contour integral in Eq. (39) can be transformed to 

Si LL \x) = -exp (£a.C F ) *(1 - x)^ sin (±a.C F t) exp (~C F e) . (41) 

The above analysis can be directly extended to next-to-leading logarithms. The jet 
function S t (x) satisfies the normalization condition Jq 1 S t (x) = 1. It vanishes at the end 
points x — > and x — > 1. The most important property of jet function S t (x) is that it 
damps faster than any power of x. 

Since the resumed factor S t (x) suppresses small x contribution, it may play crucial 
role in B decays. The B meson distribution amplitude ^(x) and the twist-3 distribution 
amplitudes 4>p{x) do not vanish at x = in general. Although the transverse momentum 
can regulate the endpoint singularity, there is still a substantial contribution coming from 
endpoint region. The factor S t (x) can suppress endpoint contribution and make pQCD 
more applicable. 

In [19], the authors discuss another resummation whose formula is similar to Sudakov 
resummation. The obtained Sudakov form factor suppresses small x contribution more 
rapidly than the factor S t (x). It shows the importance of double-log corrections from 
another point. 

The factor S t (x) presented in Eq.(41) involves one parameter integration. In order to 
simplify the numerical calculation, a simple parametrization for S t (x) is proposed [9] 
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where the parameter c is determined around 0.3. The factor S t (x) with the above simple 
parameterization form vanishes at x — 0, 1. But it does not damp faster than any power 
of x. Thus the above parametrization is proposed only for phenomenological application. 
The rigorous treatment should retain the integral. 



3 Calculations 

3.1 The formulas of B — > tt form factors in pQCD approach 

We have defined B — > ir matrix element in terms of form factors F+q in Eq. (1). The 
B — > 7r matrix element can also be expressed in another form 



^{P v )\u ltl b\B{P B )) = h{q 2 )P B , + h{q 2 )P^ 
From Eq. (1) and Eq. (43), we can obtain from / 1>2 

F? = \{h + / 2 ), Ft = \{h + / 2 ) + ^(/i - / 2 ) • 



(43) 



(44) 



7 2 



with 77 = 

In the large recoil region, the 5 — > 7r transition is dominated by the single gluon 
exchange in the lowest order as depicted in Fig. 1. The formulas for the amplitude of 
Fig. 1(a) and (b) are 



A = 



UfBm B j d £dx J d 2 l ± J d 2 k A 



x 21 ^ 7 



[XT] + l)^ B + [XT] - 1)^ B + 



(xr)m% + kj_)(£xr)m 2 B + (l± - k±) 2 ) 
m B ^x7]m B + (l±-k±y 



TTfJ, 



m B 
rriB 6 



;^b - f b )Pb„ + 2(-* B - ^b)P, 



r/ 



^B - ^b)Pb^ - 2(- - aO^Bi^ 

77 



7T/i 



- 2 



m B 6 



2 -1 + 



A-1 



xr)m 2 B + k\ ^xr\m 2 B + (l± — kj_) 2 

m B A(£;,l A 



_ fc ± ■ (Z ± - fcj + 2(Z ± - A; ± ) 2 



xrjm B + Ar[ C xr ) m B + (^l — k±) 2 I C,xr]m B + (Zj_ — Zvj_) s 



(45) 



7r/i r J 



and 



5 = 



■UfBm 2 B J d£dx J d 2 l ± J d 2 k A 



(Zrrml + P ± )(£xrim 2 B + (Z ± - k\f) 



15 




(b) 

Figure 1: Lowest order hard-scattering kernel for Bit form factor 



x 2^ T 



+ 



m B 



(46) 



m B 



Z(V b -*b) + 



£77771! + 4 ixrjm 2 B + (/ ± - A; ± ) 2 
A&Jj.),, A (l 



m B 



"(I 



k±) ■ L 



^r/m 2 B + A £xr]m 2 B + (l ± - k±) 



+ 



7] 



7T/i 



In hard scattering kernels, transverse momentum k± in the denominators are retained 
to regulate the endpoint singularity. The k\ in the numerator are power suppressed 
compared to m 2 B and they must be dropped in order to ensure the gauge invariance 
required by factorization theorem. Transform the formulas from the momentum space 
into transverse configuration 6-space and include Sudakov factors, we obtain the final 
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formulas for and Fq 77 in pQCD approach 



X 



p/tt/b^I J d£dx J b B db B b n db n a s (t) 

K) ({xt] + !)**(£, 6b) + (xv - 1)* b(£, 6b)) 



(47) 



+-^* p (x, &7r ) ( (i _ 2x)* s (£, 6 b ) + (- - 1)*b(£, 6b) ) 

%{x,K) (. 2 - \ /i^ ^ CT (a;,6 w ) 
( (1 + 2x - -)#b(£, 6b) - ^b(£, 6b) + 4 V B (£, b B ) 



-2 



m B 6 \ 
^ a (x,b v ) 
m B 6 



V 

-m B A(£,b B )) 



m B 6 



2y / C,xr]m B 



+ 



* w (x, &„) -£77(tt B (£, 6b) + tt B (f , 6b)) + 



A(C,6b)' 



+2^ p ( x , 6?r ) ((1 _ b B ) + (l + t- ^)*b(£, 6b) + 

mB V Tj m B 



hi , 



and 



pBTT 



/tt/b^b y y 6b^6 b 6 7r <i6 7r a s (i) 



(48) 



x 



* w (a:, 6^)77 ({xrj + 1)* B (£, 6b) + (xrj - l)tt B (£, 6b)) 
f^^(x, 6.) ((2 - 77 - 2x77)*B(e, 6b) + ^b(£, 6b)) 



A ^ (tf2s - 6b) - (2 - 6b)) + 4^^%^tt fl (£, 6i 

m B 6 v tbb !i 



m B 6 



-m B A(Z,b B )) 



2y/£xrjm, B 



+ 



b w ) V ^(* B (C, 6b) + * b(£, 6b)) + 



m B 



+2^%(x, K) ((n -2i + t V )* B (t, 6b) + V (l - 0*b(£, 6b) + (2 - 77)^^ 

772b V 772b 



where 



/ii = K (y/£xrj m B b B ) 



h 2 = K (y/^xrj m B b w ) 



Q(b B - h)I (y/xfj ■m B b- !r )K (^/xfi m B b B ) 
+ 9(b n - b B )I (y/xJj m B b B )KQ(^/xrj m B b^) 
9(b B - 6 7r )/ (v / ^7 ™>Bh)K (y/%rj m B b B ) 
+ 9(b n - 6b)/o(v / ^7 m B b B )K (y/£rj m B b n ) 



(49) 



(50) 
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h 3 = K_i(^/t[xri m B b B ) 9(b B - b^I^y/xfj m B b n )K Q (^/xrj m B b B ) 



+ 9(b w - b B )I (y/xrj m B b B )K (^/xfj m B b n ) . (51) 



The function Ki and 7j are modified Bessel functions with i their orders. 

The physical quantities should not depend on the choice of the scale parameter t = // 
if the calculation can be performed up to infinite orders. Therefore, the scale parameter 
can be chosen as any value in principle. However, in practice the calculation can only be 
made perturbatively. To make the perturbative expansion meaningful, the scale parame- 
ter should be chosen in such a way that can make the higher order corrections small. The 
natural choice is t — y/^xrj m B in the standard approach. If £, x — > 0, a s (t) will be diver- 
gent at t < Aq CD . When including the transverse degree of freedom, a s (/i) In y/^xrj m B j 
a s (/i) ln&e/z and a s (fj) lnfe^/i will appear in higher order corrections. Therefore we take 
t = max(y / £ai/m£, l/b B , l/b n ). The scale t = /x must be larger than Aqcd thus avoids 
the divergence of coupling constant. 

The wave functions include Sudakov corrections coming from Sudakov and threshold 
resummations 



The similar expressions are given for \l>p, ty a , \E^, fy B . The jet function S t comes from 
threshold resummation. The simplified parametrization form in Eq.(42) is taken to es- 
timate threshold resummation effects in this paper. The complete Sudakov factors S n ,B 
are given in Eq.(37). E^s are intrinsic transverse momentum dependence of pion and B 
meson wave functions. 

We compare our formulas with the results in [3, 9]. In [3], only the leading twist of 
pion is discussed. Set the twist-3 terms to zero, the two formulas of ours and [3] are the 
same except the definition of h±. The difference comes from the Fourier transform of hard 
part. In [9], the single B meson wave function ^/ B is assumed and the terms of ^/ B and 
A are neglected. The twist-3 power correction is included. The momentum projector 
in [9] for pion meson is slightly different from our projector in Eq.(7). Except for these 
differences, the formulas in [9] are consistent with ours. 

3.2 The assumption of single distribution amplitude 

Before we perform analysis of B — > n form factors in pQCD approach, we discuss an 
assumption of using single B meson distribution amplitude at first. To our knowledge, 




(52) 



S t (0exp(S B )M0Z B (b B ). 
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the use of B meson distribution amplitude in B decays firstly appeared in [20]. The authors 
suggest the simplest momentum projection for B meson which contains two terms 

Mb = -^f\i> B + m B g(0]l 5 M0- (53) 

The function <?(£) is assumed to be at the order of 1. Setting g(£) = 1, the B meson 
momentum projection reduces to 

M B = -^ B + m B } l5 M0- (54) 

This is the widely used formula in previous pQCD analysis which contains single B meson 
distribution amplitude. Here, we discuss the most familiar model used in the previous 
pQCD analysis 

<M0 = N B e(l - 2 exp (-^f ) ' (55) 

where N B is the renormalization constant makes Jq 4> B (^,)d^ = 1. In this model <\> B (£) has 
a peak at £ = A/m B where A = m B — m^. This model of B meson distribution amplitude 
has been used to fit different channels of B decays. 

B meson distribution amplitudes are important ingredients in pQCD approach. We 
should be careful about the choice of B meson distribution amplitude. In HQET, the 
definition of B meson distribution amplitude contains two terms <f> B , <p B . Compare Eq.(55) 
and Eq.(29), the difference lies in l/m B effect and can be neglected. The choice of single 
distribution amplitude amounts to taking (fi B = (fim and neglects (pm- So, the validity of 
the assumption of single distribution amplitude depends on whether the contribution of 
4>m is small or not. To test this assumption, we make approximation that <j) B ~ (j) B i = (fig 
where 0^ is given in Eq.(29). This approximation is reasonable because the calculated 
form factor F B7T (0) is: 0.239 using <p B in Eq.(55) and 0.227 using 0^ in Eq.(29). 

Table 1 shows the numerical result for the contributions from 4> B i,(f) B 2 and A B . If 
considering only <p B1 contribution, F Bw (0) is 0.227. This result is consistent with the one 
using QCD sum rule within the theoretical errors. The (fim contribution is dominant. 
For <p B 2 term, its contribution vanish at twist-3 but cannot be neglected at twist-2. The 
power suppressed term of A B is about 20% in the total numerical result. From a general 
point, there should have two leading B meson distribution amplitudes (or generally wave 
functions) in pQCD framework. But the single B meson distribution amplitude (f> B i is 
the dominant contribution. It should be noted that (p B2 contribution vanishes in the hard 
spectator scattering and weak annihilation diagrams in B — > decays. 
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Table 1: The B -> tt form factors F 5 *" = Ff5(0) with 0B2 and A B . The 
column represents contribution only. The meanings of columns 
"<t>B2" an d "^_b" are similar. The column "sum" represents the summation of 
all these contributions. 





<t>Bl 


4>B2 


A B 


sum 


twist-2 


0.042 


0.016 


0.014 


0.072 


twist-3 


0.185 





0.048 


0.233 


total 


0.227 


0.016 


0.062 


0.305 



3.3 The B — > n form factors for different models of distribution 
amplitudes 

The distribution amplitudes for pion depends on the renormalization scale. This depen- 
dence is controlled by the evolution equation. For the scale related to our discussion, 
evolution effect should be important but precise estimate of this effect depends on the un- 
known input parameters. In this paper, what we concern most is B meson wave functions 
and the reliability of pQCD framework. We will not consider the evolution effects of pion 
distribution amplitudes. Thus, the pion distribution amplitudes for both the twist-2 and 
twist-3 are taken as their asymptotic form for simplicity, i.e., 

0,,- = 6xx, (fi p = 1, (pa = Qxx. (56) 

For B meson distribution amplitudes (wave functions), three models exist in litera- 
tures: 



• Model I 



^=4--^-^ ^ (0= n-t exp( --w ) - (57) 



Model II 



<t> B \Q = — exp( ), 05(0 = — — ex P( )' ( 5 * 

where too = |A and A = m B — m h in this model. 



• Model III 



n(t,k ± ) = ^e(2t-05(k 2 ± -m 2 B m-0) , (59) 



with£ = i-. 
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Model I is proposed in [3]. It uses the equations of motion for light spectator quark 
in HQET. The ujg is at the order of Aqcd- The possible range is 0.2GeV — 0.5GeV. 
Model II is based on a QCD Sum rule inspired analysis [14]. The above two models are 
both Gaussian type and there is a peak near A/itlb- Model III uses the equations of 
motion for both light spectator quark and heavy b quark in HQET [15]. The distribution 
amplitudes have a cutoff at £ = 2£ because it takes the approximation itlb = oo 2 . The 
distributions are linear functions of £. For model I, the transverse momentum dependent 
function Eb(&b) is given in Eq.(31). For model II, the transverse momentum dependence 
part is unknown, we take the same transverse momentum dependence function Eb(&b) as 
model I. 

The other input parameters are as follows: decay constants for pion and B meson 
U = 0.13GeV, f B = 0.19GeV; A QC d = 0.25GeV; A = m B - m b = 0.5GeV; the oscillator 
parameter in the transverse momentum distribution of pion wave function (3 2 = 4GeV~ 2 
[10]; parameter in B meson wave function uj b = 0.35GeV; pion twist-3 coefficient ji n = 
2.2GeV [8]. 

We present the result of form factors F+q at large recoil q 2 = 0. Using B and pion 
wave functions and the input parameters presented above, predictions for the F^q(O) with 
the three models for B meson distribution amplitudes are listed in Table 2. 

Table 2: The numerical results of B — > 7r form factors F^q(0) with three 
models for B meson distribution amplitudes. 



Model 


0B 




A B 


total 


I 


0.564 


-0.321 


0.062 


0.305 


II 


0.641 


-0.338 


0.060 


0.363 


III 


0.540 


-0.335 


0.055 


0.260 



The predicted B — > n form factors F+jJ(0) in model I, II and model III are around 0.3 
which is favored by experiment and consistent with the prediction by other methods, such 
as QCD sum rule [11], BSW model [21]. The advantage of model III is that it only uses 
the QCD equations of motion. In reality, it is not a model. There is only one universal 
phenomenological parameter A in wave functions. 

At the end of this subsection, we would like to discuss the power corrections. Table 3 
shows that the chirally enhanced twist-3 contribution is numerically larger than leading 
twist contribution. The large twist-3 contribution is not consistent with the assumption 
2 We thank J. Kodaira for the discussions about this topic. 
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of twist expansion. The study of power corrections should be studied in a more careful 
and more systematic way which is beyond the scope of this paper. 

Table 3: The twist-2 and twist-3 contributions of pion in B — > n form factor 
F Bw = Ff S(O) with different tu B . 



u B (GeV ) 


0.25 


0.3 


0.35 


0.4 


0.45 


0.5 


twist-2 


0.086 


0.078 


0.072 


0.066 


0.061 


0.056 


twist-3 


0.345 


0.280 


0.233 


0.198 


0.171 


0.150 


total 


0.431 


0.358 


0.305 


0.264 


0.232 


0.206 



3.4 The reliability of pQCD approach in B — > n form factors 

Now, we examine the reliability of pQCD analysis in B — > 7r form factors. We choose 
model I of the B meson distribution amplitudes for illustration. The conclusions of model 
II and III are similar. 

The basic idea of pQCD approach is to use Sudakov effects to suppress the long- 
distance contribution with large transverse separations. A reliable pQCD analysis of 
B — > 7r form factors should satisfy that most of the result comes from the region where 
impact parameters b n ,b B are both small. In order to study the impact parameter b 
dependence of B — > n form facors, we introduce a cut-off b c in impact parameters b n and b B 
in the integrals of Eq. (47) and (48) by / 6jJ db B Jo* db n . In Fig. 2, we show the dependence 
of B — > 7r form factors F B7T = F B q(0) on b%. and b%. Fig. 2(a) plots the b% dependence 
with lob = 0.35GeV and fx n varied. Fig. 2(b) plots the b^ dependence with fx w = 2.2GeV 
and ujb varied. One can see that varying these two important input parameters fi n 
and uo B does not change the behavior of the F Bw dependence on b c . Similarly the b c B 
dependence of F Bn is depicted in Fig. 2(c) and (d). From Fig. 2, 50% of F Bw comes from 
the b% < 1.5GeV -1 for impact parameter b w and b c B < 2.0GeV _1 for impact parameter 
b B . The contributions from regions with large impact parameters b^.bs > 2GeV _1 are 
substantial: 31% for impact parameter b w and 50% for b B - The calculation of pQCD 
approach may still include large long-distance contributions. Therefore, the reliability of 
leading order calculations in a s expansion should be checked carefully. The more direct 
criteria is to check the distribution of the coupling constant a s (t). 

The standard to judge the reliability of perturbative analysis is that most of the 
contribution comes from the region where the coupling constant a s (t) is small. Fig. 3 
plots the form factors F Blv coming from the region where a\ < a s (t)/ii < a 2 . The last bar 
is the contribution of a s (t)/-n > 0.9. In our calculation, 70% of the result comes from the 
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Figure 2: b n and b B dependence of F Bw . 



region where ot s {t)/it < 0.2 and the contribution for t > lGeV is 38%. If we consider the 
energy lGeV is the point that the perturbation theory is broken, the non-perturbative 
contribution is about 60%. If we consider a weaker criterion that a s (t)/ir = 0.2 is the 
broken point, one can see the non-perturbative contribution is 30%. No matter which 
criterion is chosen, the non-perturbative contribution is comparable to the perturbative 
part. So, the prediction of B — > 7r form factor in pQCD approach can not be precise. 
The non-perturbative contribution constitutes the intrinsic systematic error in pQCD 
approach. 

One may also wonder if the integration over bs may not converge because the form 
factor F+ n (0), as shown in Fig. 2(c) and (d), does not saturate at any value of b c B . There- 
fore, it is necessary to study the property of the dependence of F+ n (0) on 6# in the region 
bs > 1/Aqcd- However, since the complete factor e" 5s in B meson wave function is 
divergent at Aqcd, we drop this term and perform the calculations again. We find that 
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Figure 3: The B — > 7r form factors F Bl vs the coupling constant a s (t) 

the form factors saturate at b B ~ 1/Aq C d as shown in Fig. 4. That is to say, even without 
any Sudakov suppression, the region of large bs does not contribute. 

3.5 Comparison with other works 

Whether B — > 7r form factors are hard dominant or soft dominant is a highly controversial 
topic in B physics. There are many theoretical investigations using various approaches. In 
this paper, we have studied the perturbative method. The comparison of our result with 
other recent works in pQCD approach is helpful to interpret the perturbative method. 

In [9], the authors calculate the B — > 7r form factors in pQCD approach which includes 
pion twist-3 power corrections. The intrinsic transverse momentum dependence of B 
meson wave function is considered and of pion is neglected. Neglecting the intrinsic 
transverse momentum effect in pion is unjustified. The physical reason is that Sudakov 
suppression is not so strong for realistic B decays and the large b contribution is very 
sensitive to the non-perturbative dynamics in hadrons. We have discussed that the model 
of B meson wave function chosen in [9] is nearly same as the 0j in model I. So, we find an 
explanation of this model in HQET. The reliability of B — > tt form factor in pQCD is given 
in [2] by Y. Keum, H. Li and A.I. Sanda. They observed that 97% of the contribution 
comes from the region with a s (t)/ir < 0.3. This criterion of a s (t)/7i < 0.3 is not strong 
enough to guarantee the reliability of leading order result in a s expansion series. 

The authors in [3] address some critical questions about Sudakov effects in B — > rr form 
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Figure 4: The property of F vs. b c B in B meson. The dotted curve is the result 
with Sudakov and evolution effects in B meson, the dashed one is for the case without 
Sudakov effect, i.e., without the contribution of s(£,&B,ms) in Sb, and the solid curve is 
the result without both Sudakov and evolution effects in B meson, i.e., without the total 
contribution of Sb- 

factors. We will make our comments on their questions below. 

The first question concerns theoretical issue on Sudakov form factor. For pion meson, 
it is pointed out that Sudakov form factor with next-to-leading-log (NLL) approximation 
is gauge dependent. This dependence can be cancelled by a soft function U(b) which 
resummed the single soft logarithms [22]. A complete analysis of B — > n form factors 
should include the next-to-leading-order (NLO) calculation in order to ensure gauge in- 
variance. In this paper, we take the NLL Sudakov form factor and hard kernel to leading 
order, so the gauge dependence exists in our analysis. Whether this gauge dependence 
will lead to large numerical difference is not clear at present. The NLO calculation in 
pQCD approach must include the transverse degrees of freedom, therefore the calculation 
is difficult in technique. 

The second question is about numerical dependence on the uncertainties in meson's 
distribution amplitude. This is an important and inevitable problem in perturbative 
method. We find that the three models of B meson wave functions which satisfy the QCD 
equations of motion give the consistent numerical results. More importantly, the B meson 
wave functions derived in [15] is model- independent. So, the theoretical errors caused by 
the B meson wave functions can be reduced largely. 

The third question is about the reliability of pQCD approach in B — > n form factors. 
The conclusion given in [3] is that Sudakov suppression is so weak that pQCD approach 



25 



cannot be applied in B physics. After including the intrinsic transverse momentum effects 
and threshold resummation, the perturbation behavior of the numerical result in [3] is 
similar to ours. Whether the perturbative calculation is reliable or not should be checked 
by higher order calculations in a s expansion series. Our present analysis shows that one 
can not guarantee higher order terms will be small. More works are needed before a 
completely reliable method is set up to control the higher order contributions. 

4 Conclusions and discussions 

In this paper, we have presented a systematic analysis of B — > n form factors in pQCD 
approach. The intrinsic transverse momentum effects and threshold resummation are 
important ingredients in the pQCD framework in B decays. One important finding in this 
paper is that all the three models of B meson wave functions predict F Bn about 0.3 in 
pQCD approach. This result should not be accidental. We find that the numerical results 
are not so sensitive to the choice of B meson wave functions in HQET. The B meson wave 
functions which satisfy the QCD equations of motion can reduce the theoretical errors 
largely. The consistent numerical results give us confidence that our knowledge about the 
B meson wave functions is not so poor as expected before. HQET provides a reasonable 
framework to interpret the intrinsic dynamics in B meson. 

Although the prediction of B — > 7r form factors in pQCD approach seems favorable 
in experiment, we cannot avoid many serious conceptual problems of pQCD approach. 
The calculation in our paper is at the leading order (LO). There is no strong evidence 
that NLO result is small especially in our case where most of contribution comes from 
the momentum region of t < lGeV. A critical problem in NLO calculation is how to 
deal with the dependence on the choice of renormalization scale fi. In this paper, the 
scale n is taken as the largest momentum carried by the exchanged gluon. Whether this 
choice is best or not depends on the NLO calculation. However, as we have discussed, the 
NLO calculation in B — > ir form factors is technically difficult. Another problem comes 
from power corrections. We only discuss the chirally enhanced twist-3 correction. There 
are many other power corrections, such as twist-4, higher Fock states etc. Although one 
can argue that they are power suppressed in the heavy quark limit, we don't know their 
numerical effects in the realistic case. In principle, all the above issues should be treated 
in a more systematic way than the present form. The seeming correspondence of the 
prediction of pQCD approach in B — > 7r form factors with other theoretical expectations 
is based on many assumptions and one must be careful to obtain a conclusion. To some 
extent, we think it is difficult to estimate the theoretical errors in the present pQCD 



26 



approach in B decays. 

Although there are many crucial problems in pQCD approach, completely rejecting 
it should be with caution. There are at least 30% hard contribution coming from the 
moemtnum t > lGeV in B — > 7r form factors. This contribution cannot be calculated 
reliably in any other present theoretical approaches. We cannot throw one method before 
we get the final answers. 

There is no universal, satisfactory criteria of QCD because we don't know how to quan- 
titatively calculate the non-perturbative dynamics especially in exclusive, non-leptonic 
B-meson decays at present. The complicated strong interactions in exclusive B decays 
where both perturbaive and non-perturbative dynamics contribute is a challenging task 
for the present non-perturbative methods, such as QCD sum rules, lattice QCD etc. 

In conclusion, B — > 7r form factors provide an interesting places to study the pertur- 
bative and non-perturbative strong interactions. The success of B meson wave functions 
in HQET is helpful to understand the exclusive B decays. The perturbative method pro- 
vides a simple phenomenological method to estimate the perturbative aspects of strong 
dynamics in B decays. 
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Appendix 

In this appendix we will derive the projector for B meson in the momentum space with 
including transverse momentum. We begin with the generalized Lorentz decomposition 
of the light-cone matrix element [14, 3] 



M(zU = (0\q/3(z)b a (0)\B(p B )) 



2*5(*V) + 



*- B (z\t)-n(z 2 ,t) 



75 



(60) 



a/3 
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with v = ^ and t = v ■ z. In light-cone coordinate z = (^=, z±), with z± = z ± z 3 . 



V2> y/2> 



To obtain the projector in the momentum space we consider the amplitude of one process, 

d 4 z 



T 



I 



(27T) 



-Mz)A(z), 



(61) 



where A(z) is the hard scattering kernel of the relevant process. Then 

d 4 z 



T 



M(z) I d 4 le' llz A(l) 



(27T) 4 

1 f dz + dz^d 2 z± 

(27T) 4 



M(z) J dl + dLd 2 l ± e- l( ' +Z + ^ z+ ~ r ±-^)A(l + ,LJ ± ). (62) 



For the case that A(l) is independent of Z_, we have A(l) = A(l + , l±). Then the integra- 
tion over /„ in the above equation can be accomplished to get a delta function 6(z+/2). 
Therefore we get 



T = \j d -^^M{z)\ z+=0 f dl + d*l ± e-^- r ^A(l + ,l ± ). 



(2k) 



(63) 



So we just need to consider the case z + = 0. For convenience the above equation can be 
written in the form 



T = J ^ 2 / ± i/^^e-^^-)M(,) k+=0 A(/ +) r ± ). 



Introduce Fourier transformation: 



1 f dz_d 2 z± 



(2tt)* 



e - l (^-^)^|(-il 2 , z _/2) , 



(64) 



(65) 



and substitute eq.(60) into eq.(64), we can obtain 



T 



/ 



dl + d 2 l 



-ifr 



29&l + ,h) + -f 



1 f dz_d 2 z± 

e 



it's + m B 

(2tt)= 

$iK 2 ^-/2)-^K 2 ,2-/2) 
^_/2 



^75 



A(Z+,Z ± ). (66) 



By defining A(Z + , Zj_) = J ' + dl($ B (l, l±) — ^(Z, ZjJ) and making the integrations by part, 
we can finally get 



dl+d 1 1 



_d_ 



75 



A(l+,l ± ). (67) 



To deal with we need to re-express Z M as 



28 



therefore we can get 



Substitute eq.(68) into eq.(67) and do the integration by part and drop the surface term, 
we get 



h±pL L+ {l+i l±)i+ + l±)i _ _ A(z+> w J_| 75 



•A(Z + ,Zj.). (69) 

So the projector for B meson in the momentum space is 
ifB 



M a(i = 



,I>L; + 1,1 B + *i(J+, - A(Z + , W 1 75 
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